Interaction correction to the conductance of a ballistic conductor 
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In disordered metals, electron-electron interactions are the origin of a small correction to the con- 
ductivity, the "Altshuler-Aronov correction". Here we investigate the Altshuler-Aronov correction 
SGaa of a conductor in which the electron motion is ballistic and chaotic. We consider the case 
of a double quantum dot, which is the simplest example of a ballistic conductor in which SGaa is 
nonzero. The fact that the electron motion is ballistic leads to an exponential suppression of SGaa 
if the Ehrenfest time is larger than the mean dwell time td or the inverse temperature h/T. 
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There are two quantum corrections of comparable 
magnitude to the conductivity of a disordered normal 
metal at low temperatures: the weak localization correc- 
tion, which has its origin in the constructive interference 
of electrons traveling along time-reversed paths [l], 
and the "Altshuler-Aronov correction", which is caused 
by electron-electron interactions These two correc- 
tions can be distinguished through their different depen- 
dences on temperature and magnetic field. For meso- 
scopic conductors, which are characterized by their con- 
ductance, not their conductivity, weak localization and 
the Altshuler-Aronov correction represent small changes 
to the mean conductance after taking an average over 
different disorder configurations 

Although weak localization was discovered in the con- 
text of disordered metals, it also occurs if the elec- 
tron motion is ballistic and the only source of scatter- 
ing is specular reflection off sample boundaries or ar- 
tificial macroscopic scattering sites [1, [3|- There is 
an important difference between the two cases, however: 
For weak localization in ballistic conductors a key role 
is played by the Ehrenfest time te ^S], a time which 
has no counterpart in disordered metals. The Ehren- 
fest time separates regimes of classical-deterministic and 
quantum-probabilistic motion [^, [l^, thus serving as a 
short-time "threshold" for quantum corrections in a bal- 
listic conductor. The weak localization correction ^Gwl 
has an exponential dependence on te if te is large, 
<5GvvL oc exp(— te/td — te/t^), wher e Tp is the dwell 
time and the dephasing time d, [ll|, 

In this letter we show that te also serves as a short- 
time threshold for the Altshuler-Aronov correction 6G aa 
in a ballistic conductor, in a manner quite similar to the 
way it appears in the theory of weak localization and 
other quantum corrections that do not rely on electron- 
electron interactions. In particular, we show that SGaa 
has an exponential dependence on te if te is large. 



FIG. 1: (Color online) Schematic drawing of a double quan- 
tum dot (left), with a generic set of four trajectories that 
contributes to the interaction correction to the conductance 
(right). The wiggly line represents the interaction propagator. 



to temperature is special to the interaction correction; 
Although weak localization depends on T implicitly via 
the temperature-dependence of the dephasing time, the 
corresponding exponential dependence is much weaker 
since > h/T 4\. 

For most studies of quantum corrections in ballistic 
conductors, the geometry of choice is a ballistic cavity or 
"quantum dot" . This geometry is not suitable for a the- 
ory of the Altshuler-Aronov correction, however, because 
SGaa = in a ballistic quantum dot 3, 14 1. Therefore, 
we here calculate SGaa for a "double quantum dot" , a 
device consisting of two quantum dots coupled by a bal- 
listic contact. Although SGaa will be quantitatively dif- 
ferent in other geometries, we believe that the qualitative 
features mentioned above carry over to the general case. 

A schematic drawing of a double quantum dot is shown 
in Fig. [TJ The two quantum dots are connected to source 
and drain reservoirs via ballistic contacts of conductances 
Gi, G2 ^ e'^/h respectively; They are connected with 
each other via a ballistic contact with conductance Gc- 
The classical electron dynamics in each quantum dot is 
ballistic and chaotic, with Lyapunov exponent A. (We 
assume the same Lyapunov exponent for both quantum 
dots.) The Ehrenfest time then reads 



SGaa oc g-^E/rD-a^TrE/fi^ 



(1) 



Te = A ^ \n{kpL), 



(2) 



where T is the temperature. The exponential sensitivity 



where kp is the Fermi wavenumber and L the dot size. 
We take the semiclassical limit kpL ^ 1, so that the log- 
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arithm in Eq. ([2]) is large. For a ballistic double quantum 
dot, the disorder average is replaced by an average over 
the Fermi energy or over variations of the dot's shape. 

In the semiclassical limit the electron-electron interac- 
tion in the double quantum dot takes a particularly sim- 
ple form: It is determined by the capacitances Ci,2 of the 
two dots and their mutual capacitance Cc only. (For the 
case of a single quantum dot the corresponding interac- 
tion is known as the "universal interaction Hamiltonian" 
[ist-) As a result, the interaction propagator I?(ri, r2; cj) 
for the double dot is spatially homogeneous inside each 
dot, so that it may be represented by a 2 x 2 matrix 
where the matrix indices refer to the two quantum dots. 



(3) 



with Cmn Cm^mn ^" C*c ( 1) ^ , Gmn GriiSmn -t" 

G'c(— 1)™^", and i>mn = Vm^mm where vi and V2 are the 
level densities in each dot. For all frequencies of interest 
one may neglect the first term in Eq. ([3]) and approximate 



The Altshuler-Aronov correction 5Gaa can be cal- 
culated from the interaction correction to the single- 
electron Green function Q{r,r';uj) (before the ensemble 
average). Without interactions, the relation between Q 
and the conductance G is given by the Kubo formula, 



G^^v^v^' j dydy' J d£, 



g^(r,r^;g)gA(i.^^r;0 
4rcosh2(^/2r) 



(4) 



where x and x' {y and y') are longitudinal (transverse) 
coordinates in the source and drain contacts, respectively, 
and 



VxG{-r, •)^(-,i") 



fie 
2mi 



P,^(r,-))e(-,r) 



The Altshuler-Aronov correction 5G aa then follows from 
Eq. ([4]) if one makes the substitution C/^-^ G^'^ + 
Sg^^^, with |3|,E1 



(5e«(r,r';0 = / — tanh 



dridr2a«(r,ri;0e^(r2,r';0 



X {2?R(ri, r2; a;)e^(ri, r2; e - ^) - I?^(ri, r2; ^)a« (ri, r2; ^ ^ cu)} 

I 



(5) 



and a similar expression for Sg^. Notice that Eq. ^ 
represents the Fock contribution to Sg only. The Hartree 
contribution to Sg vanishes because of the special form 
of the interaction ^(ri, r2; oj), see Eq. (|3]) above. 

Together, Eqs. dS])-© express SGaa as an integral 
over a product of four single electron Green functions. 
What remains to be done is to calculate the ensemble 
average. Hereto, we use a semiclassical approach in- 
spired by the semiclassical calculation of the weak local- 
ization correction to the conductance in Ref. [13] ■ First, 
we use the standard expression of the Green function 
Cy(r,r';^) as a sum over classical trajectories a connect- 
ing the points r and r' 18], 



^«(r,r';C) 



a^(r',r;C)* 
2tt 



(27rz?i)3/2 



(6) 



where Aa is the stability amplitude of a and Sa its clas- 
sical action. Actions at different energies are related via 



SaiO - Sai^ - Uj) = UJTa 



(7) 



where Tq is travel time from r' to r along a. Substitution 
of Eq. ^ into Eqs. (g]) and ([5]) expresses 5Gak as a sum 
over four classical trajectories ai, a2, aa, and (3 jFig. 



[U right panel]. For each set of four classical trajecto- 
ries there are four contributions to SGaa, corresponding 
to the two terms in Eq. ([5]) and the two terms in the 
corresponding expression for Sg^. The four cases differ 
by the assignment which of the trajectories ai, 012, as, 
and P correspond to retarded and which correspond to 
advanced Green functions. Each configuration of trajec- 
tories contributes to the average conductance (G) only if 
the total action difference AS is of order h systematically. 
This occurs only if the 'retarded' and 'advanced' trajec- 
tories are piecewise equal, up to classical phase space 
distances of order h^^"^ and below 0, 17]. Five possible 



configurations of classical trajectories that meet this re- 
quirement are shown in Fig. [2l In Fig. [2^, ai, a2, and 
as represent retarded Green functions, whereas (3 repre- 
sents an advanced Green function. In Fig. [2)3-e, ai and 
as represent retarded Green functions, whereas 012 and /3 
represent advanced Green functions. There are five more 
trajectory configurations that contribute to SGaa which 
can be obtained from those of Fig.[2]by interchanging the 
roles of advanced and retarded Green functions. 

In Fig.[2K, the points ri and r2 are on the trajectory f3 
(up to a quantum uncertainty), and the three trajectories 
ai, q;2, and as are equal to different successive segments 
of P (again up to quantum uncertainties). In Figs.[2j3-e, 
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FIG. 2; (Color online) Schematic drawing of five configura- 
tions of classical trajectories contributing to SGaa- Solid 
(dashed) trajectories represent retarded (advanced) Green 
functions. Five more trajectory configurations are obtained 
by interchanging the roles of retarded and advanced Green 
functions. The wiggly lines represent the advanced interac- 
tion propagator. 



the situation is more complicated because the two points 
ri and r2 need not be on /3 for this contribution to SGaa- 
If both ri and r2 are not on /?, as in Fig.[2h, the trajecto- 



ries /3 and a2 undergo a "small angle encounter" . Small 
angle encounters play a crucial role in the theory of weak 
localization and the shot noise power in ballistic conduc- 
tors d, [13, [11,0, El, • The condition that the action 
difference AS = —Sa2 +^0^ —Sp be of order fi trans- 
lates to the condition that the duration of the encounter 
be equal to the Ehrenfest time te [1, 17 1. The same is 
true if one of the points ri and T2 is on (3 and one is not, 
as in Fig. ^ and d. If both points ri and T2 are on /3, 
as in Fig. [2^, the travel time between ri and T2 must be 
less than te (20j . 

The configurations of Figs. [2^ and e are essentially 
one dimensional. Their contributions are found to can- 
cel precisely if the travel time Ta^ between ri and V2 is 
less than te. The configurations of Figs. [2b-d require a 
summation over trajectories involved in a small angle en- 
counter, taking into account the proper action difference 
AiS. For the present case, where the interaction propa- 
gator P(ri,r2;w) is independent of the precise location 
of the coordinates ri and T2 inside each quantum dot, 
this summation is essentially identical to that needed for 
the calculation of the shot noise power. Following the 
method of Refs. 0, [111, [H, one then finds 



6GaA = \ / -^H^mY.'^tni^Wi / diP^p^ 
•' m.n Wte 



H —iLot I h 



E 

k,l 



dt2dh[dr^ + 6ih) + Sit2)]P^,it2)PtPkiirE)P^'Pin{ti)e'^''''^-+'^+'^^/^ + C.C., (8) 



where T{x) = (sinhx — a::)/(cosha; — 1), P^ and P^ ~ 1 — P^ arc the classical probabilities that an electron in dot 
n exits the system through the left or right contacts, respectively, and Pmn{t) is the classical probability that an 
electron in dot n is in dot m after a time t. The first term in Eq. ([5]) comes from trajectory configurations of Fig. 
[2^ and e; the second term comes from Fig. [2)D-d. In matrix language, one has Pki{t) — [exp(— Gi/e^i>)]fe; , where the 
2x2 matrices G and i> were defined below Eq. ([3]). 

In the special case of a symmetric double quantum dot with Gi — G2 = Gd and i^i = = and neglecting the 
first term in Eq. dSl), Eq. dH]) simplifies to 

oGaa = , r,r< \3 1™ / "^TTT- uTT~ V' ^ 



where rD+ — e^vjGd and td- — e^vjiGd + 2Gc) are the 
characteristic dwell times of the double quantum dot. 
The indices -f (— ) refer to relaxation of (anti) symmetric 
charge configurations. The frequency integral in Eq. ([9|) 
can not be performed in closed form, except in asymp- 
totic limits. If Te <C min(TD± ,?i/T'), one recovers the 
result of random matrix theory [23j, 



<5Gaa 



2e2 



GdG\ 



^h, {Gd + 2Gc)3 TU+ 



-^ln^(lO) 



if tdT -C h, and 
SGaa = 



TD+ + TD- 



3T(Gd + 2Gc)3 TD+TB- ^^^^ 

if TdT 3> h. As soon as te becomes comparable to td 
or h/2'iTT, SGaa acquires a dependence on te, which 
becomes exponential in the limit of large Ehrenfest times. 



SGa/ 



2e2 GdG2(TD_e-"^/^°+ +TD+e-^E/"°-; 



{Gd + 2Ge 



(12) 
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if TD < TE < h/T, and 

2e2 GdG^TEe-^-^"^/'^ 



iJGaa = 



TT^l (Gd + 2Gc)3 
' (rD+ + ;i/27rT)(rD- + ;i/27rT) ' 



(13) 



if Te > ?i/T. Equations ((T2|) and ([TS]) reproduce the 
general exponential dependence of Eq. ([l]), with td = 
td+ taken to be the larger of the two characteristic dwell 
times. For a generic (non-symmetric) double quantum 
dot one finds the same dependence on te from the general 
expression of Eq. ([5]) . 

Let us now compare the interaction correction SGaa 
to the weak localization correction (5Gwl in the same 
system, for which one finds 



5G 



WL 



2e2 
h 



Y^dr^PlP^P^ir^) I dtPii{t), (14) 

k,l '' 



in the absence of a magnetic field. (With a magnetic 
field, (5GwL = 0.) For the special case of a symmetric 
double quantum dot this simplifies to 



5G 



WL 



2Gc(Go + Gc)^e 



2 „ — te/td-) 



(Go + 2Gc)3 



(15) 



Comparing with the expressions for JGaa derived above, 
we note that both (5Gwl and SGaa disappear oc 
exp(— te/td) in the limit te ^ td at zero temperature. 
[The fact that JGwl has a single exponential decay rate 
in Eq. IjlSp while SGaa has two exponential decay rates is 
an artifact of the symmetry Gi = G2.] At finite temper- 
ature, (5GwL and SGaa have a different dependence on 
Te, since SGwh has no explicit temperature dependence 
through the thermal time h/2TrT. 

The exponential dependence of SGaa on te/td and 
TsT/h [Eq. ^] is the main result of this letter. Although 
it was derived by explicit calculation of SGaa in a bal- 
listic double quantum dot, we should point out that the 
origin of the te dependence of SGaa is in the structure of 
the classical trajectories contributing to SGaa shown in 
Fig. [21 not in the details of the semiclassical calculation. 
It is because of this that we expect that our qualitative 
conclusions for the te dependence of SGaa carry over 
to other geometries, notwithstanding differences in the 
quantitative evaluation of SGaa in those cases [24 1. 
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